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ABSTRACT 


The stresses in the CTRAPRG and CTRIARG ring elements are not calculated 
for any of the dynamic solutions in the current COSMIC version of NASTRAN. 

This paper presents a DMAP alter sequence for Solution 8 and post-processing 
program, NASTPOST, to calculate these stresses. Test cases are presented 
which describe the method. The stiffness and the consistent versus concen- 
trated mass problems which have been ascribed to this element are reviewed. 

The DMAP alter sequence introduces Solution 8 displacements to a Solution 
1 module to calculate Real and Imaginary stress components during the execu- 
tion of Solution 8. The post-processor, NASTPOST, calculates the magnitude/ 
phase stress results. 

The DMAP sequence has been written specifically for Level 52 f^C/NASTRAN, 
but can certainly be used for any COSMIC version with slight modification. 


INTRODUCTION 


None of the currently documented versions of NASTRAN calculate the 
dynamic stresses in the CTRAPRG and CTRIARG solid of revolution elements. The 
stresses for these elements are calculated in NASTRAN for static solutions 
(e.g.. Solution 1) but not in the dynamic solutions (e.g.. Solution 8)._ Com- 
ments have been made by others which express the reasons for not including the 
stress calculations are related to the formulation of the mass matrix for the 
element. 

Sample problems are given to show that the difference between the consist- 
ent and concentrated mass approach is greater than one might expect from argu- 
ments solely between the merits of consistent or concentrated mass. 

This paper describes a DMAP alter sequence for Solution 8 and a post- 
processing program, NASTPOST, to calculate these dynamic stresses. The DMAP 
alter sequence introduces the displacements computed in Solution 8 to a Solu- 
tion 1 module to calculate the complex stresses in the form of real and 
imaginary components. The post-processor, NASTPOST, calculates the stresses 
in the form of magnitude/phase. 
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DISCUSSION 


It is not spelled out in the NASTRAN Users Manual that stresses for the 
solid of revolution elements are not calculated for dynamic solutions» There- 
fore, if one asks for stresses in a Solution 8 case control, the run is not 
aborted, but no stresses are obtained. 

In order to perform noise path studies of an axi symmetric structure it 
became necessary to obtain these stresses. At first, the displacements for 
the entire structure, obtained from a Solution 8 forced vibration analysis 
were written into an output file; then these displacements, less one, were 
written into SPC format as enforced displacements for a Solution static analy- 
sis (this was done for the real and imaginary components separately). This 
technique was later modified, utilizing the DMAP alter sequence A0S8$CS and a 
post-processor, NASTPOST. 

The DMAP alter sequence is given in Figure 1. The major points are: 

• The user can specify output requests as usual for SPCFORCES 
and DISPLACEMENTS. 

• The user should specify STRESS (PUNCH) = ALL or a particular 
set ID if he wishes to subsequently use NASTPOST to calculate 
the magnitude/phase. This punched file will be sent to the 
users system space. (FOR 013. DAT for the MSC/NASTRAN VAX 11/780 
VERSION). 

• A0S8$CS should be pldced on the user's RFALTER library and 
executed then by calling RFAI = A0S8$CS. 

The program NASTPOST is given in the appendix and is used to calculate 
magnitude/phase stress components from real/imaginary stress components. The 
major points are: 

• The components from F0R013.DAT above, are used as input to 
calculate the magnitude/phase stress components. 

• This program can be run immediately after the execution of 
MSC/NASTRAN or at some later time. 

The test problem for A0S8$CS and NASTPOST is a circular plate fixed at 
the edges and driven by a single force, 100 dynes, at the center, normal to 
the plane of the plate. The finite element control model is the CQUAD2 and 
CTRIAG2 bending element model shown in Figure 2. The CTRAPRG model, shown in 
Figure 3, is formulated as a concentrated or consistent mass for each of the 
runs. The NASTRAN default value is the consistent mass matrix. The concen- 
trated mass matrix is entered as C0NM2 data. The three cases are compared in 
Table 1 for static, 2000 Hz and 8000 Hz at a position near the concentrated 
load and at the fixed edge. 


64 



The concentrated mass formulation gives good results, as compared to the 
control model. The consistent mass, or default formulation, gives results 
which do not agree with the control model at either the low, 2 kHz, or high, 

8 kHz, forcing frequencies. 

The static solution agrees very well with the control model which indi- 
cates that the stiffness of the model is represented correctly by solid of 
revolution elements. The error therefore is associated with the mass matrix 
formulation. The degree of error is obviously greater than one would expect 
from the normal arguments of consistent versus concentrated mass differences.^ 

It can be argued that the use of cyclic symmetry with 3D elements rather 
than solid of revolution elements would have been a possible solution. This 
is certainly an avenue that deserves added investigation for comparison of 
cost and accuracy of solution compared to the solid of revolution elements with 
concentrated mass matrix. 


CONCLUDING REMARKS 


A DMAP alter sequence for Solution 8 and a post-processing program 
NASTPOST has been presented to calculate the dynamic stresses in CTRAPRG and 
CTRIARG solid of revolution ring finite elements. Users of this technique are 
cautioned to use the concentrated or lumped mass matrix rather than the con- 
sistent mass (default value) matrix. 

The DMAP sequence has been written specifically for Level 52 MSC/NASTRAN, 
but can certainly be used for any COSMIC version with slight modification. 
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TABLE 1 


COMPARISON OF STRESSES, 3/8 cm from CONCENTRATED LOAD 


FREQUENCY 

01 

2 kHz 

8 

QUAD2 

134.4 

75.5 

66.4 

TRAPRG (CONS.) 

132.3 

17.2 

63.1 

TRARG (CONC.) 

132.3 

96. 

60.5 


TABLE 2 

COMPARISON OF STRESSES, 3/8 cm from FIXED EDGE 


FREQUENCY 

01 

2 kHz 

8 

QUAD2 

44.4 

34.2 

38.2 

TRAPRG (CONS.) 

45,6 

27.0 

10.0 

TRAPRG (CONC.) 

45.6 

33.0 

36.0 


1 OBTAINED FROM SOLUTION 1 

6 £ 



FIGURE 1 - ALTER A0S8$CS 


BEGINNING OF ALTER A0S8SCS 

THIS ALTER PACKAGE IS USED TO CALCULATE 

^DISPLACEMENTS ( REAL/ I MAG I NARY ) OR 

(MAGNITUDE/PHASE) 

tSPCFORCES (REAL/IMAGINARY) OR 

(MAGNITUDE/PHASE) 

^STRESSES (REAL/IMAGINARY) 

FOR THE CTRAPRG AND CTRIARG RING ELEMENTS 


CASE CONTROL INPUT 



FIGURE 1 - (Conf d) 


THE USER SHOULD SELECT THE DESIRED 
OUTPUT AS USUAL FOR DISPLACEMENTS 
AND SPCFORCES. 

THE USER SHOULD SELECT THE PUNCH 
OPTION FOR STRESS IF IT IS DESIRED TO 
SUBSEQUENTLY CALCULATE (MAGNITUDE/ 
PHASE) USING A POST-PROCESSING PROGRAM 


ALTER 166 

OFP OPPCl,OQPCl,OUPUCl,,,//U^N,CARDNO t 

ALTER 185,186 

PARAM //STSR/13/-64 f 

GP3 GEOH3,EQEXIN,GEOM2/,ETT/O/U,N,NOGRAU/0 S 



FIGURE 1 - (Cont'd) 


cn 

>X) 


PARAHL UPVC//C,M,TRAILER/a/y»N»ROUS * 

MATGEN ,/UNIT/l/ROUS t 


nOOTRL UPUC////3 * 
nPVAD IWIT,UPUC,/ASQR/ « 
DIAGONAL ASQR/ATRH// * 

ADD UPUC,/BSQR/(0.O.-1*®> * 


DIAGONAL BSQR/BTRH// S 

SDR2 CASECC,CSTN,HPT,DIT,EOEXIN,SIL,ETT,EDT,BGPDT,,, 

XVCDB/» , ,OESCR» »^STATICS/S,N,N0S0RT2 * 

SDR2 CASECC,CSTn,nPT,DIT,EQEXIN,SIL,ETT.EDT,B6PDT,,, 

XVCDB/, , .OESCI „/STATICS/S,N,N0S0RT2 « 

OFP ,,,OESCR,,//S,N,CARDNO % 

OFP ,,,OESCI,»/-'S,N,CARDNO % 


ATRN,EST 

BTRN.EST 


» 


» 


PARAH //STSR/7/-64 f 


ENDALTER f 

f 
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APPENDIX A 


THE NASTPOST PROGRAM 
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c DATA SET NASTPOST LEVEL 

COlWON 'MDRCOH/TITLE(lfi),SU»T<t6).W^ 

DATA DTI T/ ' tT IT ' CASE/ 'CASE , ‘fUD /, 
t DELE/'»ELE'/,BSTR/' STR'/,DLA»/ »LAR / 
DATA 1036.1037/3*6/ 

1 CONTINUE 
REMIND 7 

C - GET TITLE CARD 

5 CONTINUE ,T«c 

REA D<7. D66.END;8a9) TENP.TITU 
IF(TW.EQ.DTIT) CO TO S 
GO TO S 

C - GET SUDTITLE CARO 

® READ<7%.g»tt;8”> IS"?'*"” 

IF(TEW.EO.DSUD) go TO 7 
GO TO C 

C - GET LABEL CARD 

’ HEADcffloO.END-DW) TEHP.LAiEL 
IF(TW.EQ.DLAD) go to 10 
GO TO 7 

C - GET STRESS CARD 
10 CONTINUE 

HEAD<7,«16.Er«-»99) TERR 
IF(TEHR.EO.DSTR) GO TO DO 
GO TO 10 

e - GET SUDCASE IDENTIFICATION 

** »S<y936.END"WG) JO^i**** 

IF(TEV.CG.CASE) GO TO 30 
GO TO 36 

C > GET ELENENT TYPE 

*• «?l55wD.wp.^ 

IF(TEV.NE.DELE) GO TO « 

C - CteCE ELENENT TWES 

IF(IELTVR.EG.M) «« « 

IF<IELTVR.EG.37» GO TO 370 
GO TO f __ _ 

C - ELENENT TYPE • 3S 

^ «?f]isy G) ORlt RU3G(XS1D,XELTYP.1E0F) 
IF(X03S .Et. II CALL 

IFlim .EG. 1 .AND. XEOF .EG. II GO TO ODD 
X03I ■ H0OXXO3BTt.ai 
GO TO G _ 

C - ElflKNT TYPE > 37 

^ T^flSnf .EG. 61 CALL RIOTtXSXD.XELTYP.XEOFI 
XF(X037 .EG. II CALL RC37( XSXD. XELTVPtXEOF I 
XFIIW .EG. 1 .AND. XEOF .EG. II GO tO ODD 
X037 - N0Da037*l,8l 
GO TO G 


900 F0NNAT(A4,SX.1SA4«AI} 
910 FCRNATCDK.AG) 

936 F0RHAT(4X.A4,9X*X9I 
936 


•} 00<93 

00004 

00005 

00006 
00007 


00009 

00616 

06611 

06613 

06613 

06614 

06615 

06616 
06617 
0661D 
06619 

06636 
06631 
00633 
06633 
66634 
00635 
66636 

06637 
0663D 
0M39 
06636 
66631 
66633 

66633 

66634 

66635 

06636 

06637 
66630 
66639 


60641 

00043 

06643 

60044 

66645 


60047 

66040 

66649 


66651 

66653 

66663 


66655 

66656 
00657 
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i<d>3S9 


F0RWAT(A4,12X.ni) 

END 

C DATA SET NASTRU36 AT LEUEL 004 AS OF 11/02/79 

SUBROUTINE RU36( ISID. lELTVP, lEOf ) 

OIHENSION TEnP(a).OATA<4J 

DATA TITUE/'fT V.CONT/' -CON '/.BLANK/" '/ 

DATA INN.IOUT/7.9/ 

REWIND lOUT 
PRINT 10 

10 FORNAT< 'SUBROUTINE RU3£') 

READ! INN, 900, END*999 ) lELNO, DATA! 1 ) , DATA < 2 ) . DATA ( 3 ) 

001 CONTINUE 

READ (INN, 910, END-990) CARDN,DATA(4) 

IF(CARDN .NE. CONT) <K> TO 990 
URITE(IOUT) ISID,IELTVP,IELNO,DATA 
C READ(INN,920,ENO-999) TENP 

C BACKSPACE INN 

CALL BACKSPI TENP, INN. 1999) 

IF(TENP(l).EO. BLANK) 

9 REAO(10,900.ENO-999)IELNO,DATAa).DATA(2).DATA(3) 
IF(TEHP(1) .EQ. BLANK) GO TO 001 
IF(TCNP(1) .HE. TITLE) GO TO 990 
800 CONTINUE 

ENOFILE lOUT 
REUINO lOUT 
RETURN 

990 CONTINUE 
STOP 3600 
999 lEOF • 1 
GO TO 800 

900 FORNAT(I10,8X.3E18.6) 

910 F0RHAT(A4,14X.3E18.«) 

920 F0RHAT(2A2) 

END 

C DATA SET NASTRW37 AT UWEL 004 AS OF 11/02/79 

SUBROUTINE RU37(IS10.1ELTVP.1E0F) 

10 FORNATC 'SUBROUTINE RU37' > 

DIHENSION TEHP<2),OATA(20).KKREAD(33) 

DATA TITLE/'9T '/,CONT/'-CON'/,BLANK/' '/ 

DATA INN,I0UT/7,8/ 

REUIND lOUT 
PRINT 10 

READ! INN. 900. END-999 ) lELNO.DATAt 1 ). DATA(2 ) , DATA( 3 ) 

001 CONTINUE 

REAO(INN.910.E)W-990) CARDN.0ATA(4).DATA<5),DATA(6) 
IFiCARON .NE. CONT) GO TO 980 

REAO(1NN.910.END>990) CARDN.0ATAC7).DATA(8),DATA(9) 
IFICARON .NE. CONT) GO TO 990 

REA9aNN.910.ENO«990) CARDN.DATA(10).DATA(11 ).DATA(12) 
IFICARON .NE. CONT) 00 TO 9M 

READ! INN. 910.END-990) CAR0N,DATA(13).DATAtl4).DATA(15) 
IFICARON .NE. CONT) 00 TO M 

READI 1NN.910.ENO-990) CARON.DATAI IS ).DATAI 17 ).DATAI 18 ) 

IFICARON .HE. CONT) GO TO 9M 

READI INN. 910. ENO-990) CARB^DATAa9).DATA(20) 

IFICARON .NE. CONT) 00 TO m 
URITEIIOUT) ISIO.IEITV^IELNO.OATA 
C READI INN. 920. ENO-999) TENP 

C BACKSPACE INN 

READ(INN.930.ENO-999)KKREAD 
REUIND 10 

URITEI 10.930 )KKREAD 
REUIND 


• -S 


dae<^2 

00004 

00005 


00007 

00008 

00009 

00010 
00011 
00012 


00013 

00014 

00015 

00016 

00017 

00018 

00019 

00020 
00021 
00022 

00023 

00024 

00025 

00026 

00001 


00003 

00004 

00005 

00007 

00008 
00009 
09010 
00011 
00012 

00013 

00014 
90015 
00016 

00017 

00018 

00019 

00020 
00021 
00022 
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!• 


REftD(l«.92®>TEnP 
REWIND !• 

t^Kplael ^!eAd-999^ i £ L.O. DATA( 1 ) . DATA ( a ) . DATA ( 3 ) 

* IF (TEHPdl.EO. BLANK) <iOTO ®®* 

IF(TEHP<1) .HE. TITLEJ '^0 TO 99* 

8M CONTINUE 

ENDFIU lOUT 
REWIND lOUT 
RETURN 
99« CONTINUE 
STOP 37RR 
999 lEOF • 1 
CO TO 8M 

9M roRNATdio.sx.^ta.fi) 

919 FOHNATCA^.HX.aEia.fi) 

989 F0RNAT(2Aa) 

939 FORMAT! 33A4) 

^ r.R0S?I5EiJ3E^n&^^^ 

PRINT 19 

lELCNT -99 

KA0<lNii.999fEND?999) lELNO.DATAK 1 ),DATAI(2).DATAI(3) 

5i;g??g5Ti;;!siMi» 

IFdELNOR .NE. lELNO) CO TO 999 

•£■«!«!<!! -S- :;5i ?54l.!?* 


DATAKDBDATAKD) 


IF(DATAKl) .EO. •••» LTSSiii 
IFlOATAId) .«T. 9.9) PHASEd) 

IMDATAId) .LT. 9.9) PHASEd) 

90 TO C99 

RATIo'^AWlOATAId I/MTABUi’ 

* IF<^I(I).t T.9.J ..yO.,PATyil).lT.9.9) 

** IF^icnll^ ?AHO. STARd>.9T.9.9) 
X PMASCd) • PHASEd) ♦ ET9.9 

e **• SlTEdPRT.939) 1SI0.IE^ 

IFdELGNT .LT. 59 ) 00 TO 799 
CALL HKM11SI9) 
lELCHT « 9 

gjsri Ku*. . . 


9.9 

99.9 

279.9 


<•0054 

O0055 

00856 

00057 

00058 

00059 

00030 

00031 

00032 

00033 

00034 

00035 

00036 

00091 

00092 
00003 
09004 

00005 

00096 

00008 

00009 

00010 

00911 

00912 

00913 
09914 

00915 

00916 
90917 

99918 

99919 
99929 
99921 
09922 
00923 
09924 
9992S 
09026 
00927 
09928 


09931 

00932 

09933 

00934 

00935 

00936 
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URITE(IPRT.94«) IELHO.‘<RHgG<I).PHASEa))>I-l.‘*> 
REA0(IN«,9a®,END*989) TEHP 

’"cSLllcitSPtTEIlP.IllH.4SM> 

. SS!::i«“HliS!S!.E^.»5T..ti..sp.«.<E>.»«'«>3> 

*IF(TEHP(1) .M. bl^) fiO TO w 
IF(TEnP(l) .N£. TITLE) CO TO 890 
RETURN 
999 CONTINUE 
STOP 3691 
999 lEOF • 1 

999 F0MMT(I19<8X<3E18.6) 

919 F055TIM.i'»X.3£18.8) 

^ :is?I!!E’SsSnELj!.tl5\ 

bStA IP»T.INN,I0UT/6,7.8/ 

PRINT 19 
lELCHT « 19 

ReSo<1N;.I99!eND? 999) IELN0.D9TRm),BATAI(2).DPTAI(3) 
REa2<INn!1i9,EN^> 

S^iS??»p9f9:ESS%|“c52^^ 

Ji;S??SS.9f9:ESWc»TAI<« 

IM1SI»\«. ISIW^gO TO 9M 
irUEtTVP .HE. lEtTPR) CO TO »9 
IFUELMW .re. lElHO) CO TO 999 

"■K<‘ 

SiSSS!!! :S: t;: M.ff; « 

S!K:!!5! :«: *.::! SSSii! • 

00 TO 090 

Sno"^i'<«S<BRT9I<I )/BRT9R<I ) » 


DATAKDTWTAia)) 


oa038 

«i9039 

08040 


00041 

00043 

00043 

00044 

00045 

00046 

00047 

00048 

00049 

00050 

00051 

00052 

00053 

00001 

00002 

00003 

00004 

00005 

00006 
00008 

00009 

00010 
00011 
00012 

00013 

00014 

00015 

00016 

00017 

00018 
00919 
00029 
00021 
00922 
00023 
00924 

00025 

00026 
09927 
00928 
00029 
09939 

09931 

09932 
00933 

00034 

00035 

00036 
09937 
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X PHASE(l) ■ PHASE(I) + 18*.» 

IF(DATAI(I).LT.a.« .AW). DATAR(1 ).GT.«.») 

X PHASE(I) • PHASE(I) * 27*. • 

699 CONTINUE 

URITEaPRT,93«) ISID.IELTYP, lELMO.DATAR.DATAI 
IFUELCNT .LE. 7) CX) TO 7M 
CALL H037(ISID> 
lELCNT • 0 
7«0 CONTINUE 

lELCNT • lELCHT ♦ 1 
DO 710 I - 1.5 
J • ■IXll-l) ♦ 1 
K • J 3 

IF(1 .EQ. 1> URITE<IPRT,940) lELNO.I. 

X ((RM)fi(IXl>.PHASE(lXl>>.IXl*J.K) 

IF(I .NE. 1) URITE<IPRT,9S0) 1, 

X ((RNAG<1X1).PHASE(IX1)),1X1-J,K) 

710 CONTINUE 

URITEdPRT.MO) 

: REA0(INN.920.EN0*9S9) TENP 

: MCKSPACE INN _ 

CALL DACKSP<TENP.INN.i999) 

« KAO( lolsM. EHD*999 ) lELNO. DATAI ( 1 ) , DATAl < 2 ) . DATAl ( 3 ) 
If(TEHP<l» .EO. BLANK) GO TO Ml 
1F(TEHP<1) .NE. TITLE) GO TO 990 
RETURN 

999 CONTINUE 
STOP 3791 
999 lEOF • 1 
RETURN 

9M F0RHAT(I19.8X,3E18.6) 

919 F0RNAT(A4. 14X.3E18.fi) 

Qfift F0RMTC2A2) 

938 F0RHAT(lX.3I18.ie</.4(SX.lPE13.fi))) 

948 F0RNAT<1X.I5.1X,I3.4X,4<1PE12.5^^'.8PF18.5.5X)) 

958 F0RNAT(7X,I3.4X,4<1PE12.5.' /'.8PF18.5.5X) ) 

968 FORNATC ' ' ) 

END 

C BATA SET NASTHD36 AT LEUEl M7 AS OF 18/24/79 

SUBROUTINE HDSfidSIO) 

18 FORNATI 'SUBROUTINE HB3fi'> _ , 

CORHON /H0RC0tl/TlTU(16).SUBTd6).LABELd6> 

PRINT 18 
IPRT ■ 8 

URITEdPRT.lM) TITLE 
URITEdPRT.118) SUBT 
URITEdPRT.lM) LABEL, ISIS 
URITCdPRT.148) 

URITE (IPRT, 158) 

URlTEdPRT,lfi8) 

URITEdPRT,178) 

RETURN 

IM F0RIIAT('1',3)(,16A4,AR) 

118 FORHAT(' ',3X,1SA4,A8) _ 

IM F0RHAT('8',3X,15A4,A2.S8)(, 'SUBCASE', 13) 

139 FOiiMITC * ^ ) 

148 F0RNAT(27X,'S TRESSES FOR THE TRIAN 


3M38 

>)d039 

00341 

00048 

00843 

00844 

08845 

08846 
03847 

08848 

08849 
08858 

08851 

08852 

88853 

88854 

88855 

08856 

08857 
0885B 


88859 

88868 

08861 

08862 

88863 

88864 

08865 

08866 

88867 

88868 
00069 
88878 
80871 
88072 
88873 

80M1 

88082 

88M3 

88884 

88885 

88886 
88M7 
88988 
88889 
8M18 
8M11 
8M12 
8M13 
8M14 
8M15 
8M16 
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X 'GULMR RINGS (CTRIAR6)') 

15« FORHATCGlX.'CnAGNlTUOE/PHASE)') 

16G F0Rf1AT(4X. 'EL'.17X. 'RADIAL', ISX.'CIRCUNFERENTIAL', 

X 19X, 'AXIAL'. 24X, 'SHEAR' } 

17# F0RHAT<4X,'ID',19X,'<X)'.24X,'<TMETA)'.24X,'(2)', 

X 26X.'(ZX)' ) 

END 

C DATA SET NASTHD37 AT LEVEL «A6 AS OF lA/24/79 

SUBROUTINE HD37(ISI0) 

IG FORMAT! 'SUBROUTINE HD37' ) 

COmON /HDRC0n/TITU<16),SUBT(16).LAB£La6) 

PRINT IG 
IPRT • 6 

URITE(IPRT.IGG) TITLE 
URITECIPRT.llG) SUBT 
URITE(IPRT.12G) LABEL. ISID 
URITE(IPRT.14G) 

WRITE! IPRT. tSG) 

URITEUPRT.tfiG) 

URlTEdPRT.lTGI 

RETURN 

IGG FORMAT! 'l'.3X.lSA4.Aa) 

IIG FORMAT!' '.3X.1SA4.A2) 

12G FORMAT! 'G'.3X.1SA4.A2.5GX. 'SUBCASE'. 13) 

139 FORHAT( ^ ^ ) 

14G F0RMAT!27X.'S TRESSES FOR THE TRAPE', 
X 'ZOIOAL RINGS (CTRAPRG)') 

ISG FORHAT!G1X.'!MAGNITUOE/PHASE)'> 

16G FORMAT ! 4X , ' EL ' , 2X. ' ST ' . 13X. ' RADI Al ' . 1 9X. ' C I RCUMFERENTl AL ' , 
X 19X. ' AXIAL'. 24X. 'SHEAR') 

17G FORMAT! 4X, ' ID' , 2X. 'PT' . 15X. ' !X ) ' . 24X. ' ( THETA ) ' , 24X. ' !2 ) ' , 
X 2EX.'!ZX)') 

END 

SUBROUTINE BACKSPiTEMP.INN.t) 

DIMENSION KKREA0!33).TEMP!2) 

READ! INN.93G.END*999)KKREAD 
REWIND IG 

WRITE! IG.93G1KKREAD 
REWIND IG 
REA0!1G.92G)TCMP 
REWIND IG 

93G F0RMAT!33M) 

DM FORNAT!2Ra) 

RETURN 

909 RETURN 1 

END 


O0G17 

OGGIB 

09019 

00020 

00021 

00022 

00023 

09001 

000G2 

00003 

00004 

00005 
0000S 
90007 
00008 
00009 
0O01G 
00011 
00012 
90013 

00014 

00015 

00016 

00017 

00018 
00019 
90020 
00021 
00022 
90023 



